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ABSTRACT 



Computation of the terms of the vorticity equation 
is made at 300 mbs using the assumption that vertical mo- 
tions are negligible and approximating relative vorticity 
by geostrophic vorticity computed from finite differences. 

The validity of the computations is then tested, giving 
good qualitative agreement but poor statistical agree- 
ment.. The isobaric divergence at 300 mbs as computed 
using is compared with the value of 

divergence as obtained from the vorticity equation. Again 
there was only qualitative agreement. The magnitude of the 
lateral and longitudinal components of isobaric divergence 
were studied and it was found that the lateral component 
was dominant in 75 per cent of the cases studied. No rela- 
tionship was found between the longitudinal component and the 
total divergence. The synoptic aspects of the divergence 
fields at 3OO mbs are discussed with respect to the 500-mb 
surface. 

The witer wishes to express his appreciation for the 
guidance and encouragement given him by Professor F. L, Martin 
of the U.S. Naval Postgraduate School in this investigation. 
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- Rectangular (Cartesian) coordinates; standard 
coordinates with x axis toward the east, y axis 
toward the north, z axis toward the zenith. 

- Rectangulat (Cartesian) system of unit vectors 
along X, y and z axes. 

- Vertical component of absolute vorticity. 

£ A three dimensional time rate of change operator. 

- Gradient operator on an isobaric surface. 

- Wind velocity vector. 

- Isobaric divergence 

- ^ - Rate of change of pressure of an air 
^ particle with time. 

--Local rate of change operator \d.th respect to 
pressure. 

- Horizontal component of wind vector 
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Time rate of change operator along an isobaric 
surface. 
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u,v,w 
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- Acceleration due to gravity (32.17 feet per second 
per second). 

- Coriolis parameter (2 fl sin jp ). 

- Angular speed of the earth (7.292 x 10“^ radians 

per second). 

* 

- Latitude ' 

- Height of an isobaric surface. 

- Laplacian of height. 

- Components of wind velocity along the x,y and z 
axes respectively. 

- Space-aver aged-height of an isobaric surface. 

- Grid distance (2.5^ of latitude -9.12 x 10^ feet) 
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r. ^ “ Rate of change of cor lolls parameter along 
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the y axis. 

- Angle between wind velocity and the gradient 
of (z-z). 




- Density of the air. 
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- Natural system of unit vectors. 
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- Wind direction. 
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- Correlation coefficient. 
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- Covariance of stochastic variables x and y. 
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1. Introduction. 



Several studies of divergence in the atmosphere have 
been carried out by different methods and by different inves- 
tigators, Some of the most notabld of these are Fleagle ^2] , 
Arnason [l] , and Landers [5^ • Their investigations have 
involved a three dimensional study of the distribution of 
divergence in the atmosphere. Landers [5] computed horizon- 
tal divergence using constant level charts and a "del com- 
puter" which consists of computing divergence from the wind 
vectors at the vertices of an equilateral triangle with an 
altitude of of latitude. Arnason [l] computed the inte- 
grated divergence over a 12 hour period in isobaric layers 
with the use of vertical velocities, the local rate of 
change of thickness and the temperature advection, Fleagle 
computed horizontal velocity divergence and mass divergence 
up to about 150 mbs and found a level of non-divergence at 
about 600 mbs. 

This study is carried out at 300 mbs only and no attempt 
is made to extend it to ttoee dimensions. The main objecti-'^e 
of this study \ra.s to gain some practical knowledge of the 
behavior of divergence patterns arjd to apply the results to 
the prognoses of the 500-mb chart. Computations of the di- 
vergence and the individual rate of change of the absolute 
vortlcity with time are made, and their relationship tested. 
Also a study of the relationship between the longitudinal 
component of the divergence and the isobaric divergence is 
made . 
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2. Source of Data. 



A series of 300 -mb maps, carefully analyzed by the 
staff of the U.S. Naval Postgraduate School, was used in 
making the computations. The analyses employed were those 
of 150OZ, 1 ^ November, 03OOZ and 150OZ, 16 November 1955 * 
On the first chart there is a trough and a closed low 
at about 112 ° W and a ridge at about 84 ° W \^th consider- 
able meridional flow. At 0300Z on the l6th the trough and 
ridge have moved about ^ ea?tward, retaining about the 
same intensity... By 150OZ on the I 6 th.the amplitude of the 
trough-ridge system has decreased considerably and moved 
to 98° W and 72 ° W longitude so that in between the trough 
and ridge the flow is becoming zonal.. Computations, which 
will be described in the following sections, were carried 
out betv/een this trough and ridge. In addition some com- 
putations were made behind the trough on the 0300Z, 16 
November chart,. A total of I38 computations of all five 
terms of Equation ( 7 ) below were made, 

3. Mathematical Analysis, 

The vorticlty equation in isobaric coordinates is: 

( 1 ) 4r = - is- V'pw 4^' 

dZ ~ a P . 

This is not in the form in which computations from a 
synoptic chart can be readily madei Since vertical veloc- 
ities are normally small, the second term on the right side 
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of the equation may be ignored.. Equation (1) then becomes 
the approximation to the vorticity equation as it is common- 
ly made in meteorology. 

( 2 ) 

Making the geostrophic approximation, J. SJ^z ? and 
noting that the coriolis parameter does not change v/ith 
time, Equation (2) becomes: 



In order to compute the individual rate of change of 
vorticity it would be necessary to construct isobaric tro- 
jectories of the air parcels. However, such computations 
may be avoided by making use of the fact that; 



( 4 ) 



Dt 



Air 



+ y„ • ypi 



Equation (3) then becomes; 

Since £ is qua si- constant compared to the horizontal var- 
f 

iation of \7 z , it may be brought outside of the differential 
operator in the second term of the left side of Equation (5)« 
In order to get Vz in a fomn in which computations can be 
made, finite differences nay be used to give; 



(6) V’-z = ^ 



3 



Combining Equations ( 5 ) and ( 6 ) and multiplying by fd ^ gives: 




^ d e 



The equation is now in a form in v/hich computations were 
made. The left hand side of this equation now represents 
fd^D ia, and the right hand side represents; -fd^j’^f V«' V). 

By measuring wind velocities in feet per second and grid 
distances- in feet, the units of this equation can be shown 
to be feet per second. 

Computational Procedure. 

The quantity z may be readily computed using a grid 
as shown in Figure 1. 

U 

Figure 1 

Grid for computing space mean. 

The value of 2d was chosen to be 5° of latitude. This 
value ix about half of the length of the short-wave trough 
so that in constructing a "Z-chart, the short-wave troughs 
tend to be removed, but the long-wave troughs are xmaf- 
fected. 
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The construction of the z-chart where z is’ z^t Z2+ 23+ zl^. 



(Ji'igure 1) , was accomplished by displacing two z-charts 10 
degrees of latitude laterally and adding graphically giving 
z^+ Z3. The t\ro maps were then displaced 10 degrees of 
latitude northward and added giving Z24z^.. The 23^4 Z3 
contours are then added to the Z24zi^. contours giving h'z. 

By dividing by 4 we then have a space mean map from which 
the actual map may be graphically subtracted giving iso- 
lines- of (“Z-z). The local rate of change of (z-z) was 
approximated by the average rate of change in the quantity 
over a 12-hour interval, as measured at any given station. 

The advection term [term .(b^ was computed by measur- 
ing (1) the gradient of (z-z) across selected stations 
over a distance of 5° of latitude, (2) the vrind velocity 
in feet per second, and (3) the angle between the -wind 
vector and the gradient vector.. The magnitude of term (b) 
is then |Vl (Vp(z-z)l \cos i/^| . a diagram of the vectors 

involved is shown in Figure 2. 




V,(2-z) 



\J 



Figure 2 



> 



Diagram of Angle ^ 



AJLl factors of term (c) with the exception of v are 
constant for a given latitude and can be computed for each 

5 



latitude circle. Taking the appropriate value of fd^ ^ 

g 

and multiplying hjr the northward component of the wind 
velocity gives the final term in the vorticity computation. 
Now by equation (7) the si^n of terms (a)-, (b) and (c) rep- 
resent the quantity fd^ D Ja, subject to the accuracy of the 

g Dt 

approximations made in neglecting vertical motions and in 
computing (z-z). 

Several methods have been devised for computing iso- 
baric divergencd [b> 6^ • It was decided to use the fol- 
lowing equation in the computation. 



( 8 ) 



^ u + Ay. 



A square grid 2^ degrees of latitude on a side was used. 

An average u was computed on each side of the square and 

an average v at the top and bottom. As shown in Figure 3> 

the isobaric divergence is then U 2 -U]_ j_V 2 -V 2 ^ 

d ^ d 





Figure 3 

Grid for computing divergence. 



Since the grid is square no spherical correction is 

necessary, and grids can be oriented along the major flow. 

If the values of isobaric divergence are then multiplied 

6 



by tiarm (d) of Equation (7) the value, should be the same 
as the sum of terms (a), (b) and (c). The divergence was 
also computed by dividing the sum cxf terms (a), (b) and (c) 
by tenn (d) of Equation (7)«- For the discussion in this 
paper this value of isobaric divergence will be referred 
to as effective divergence . 

Isolines of the computed values of fd ^ D 'S<y are 

g Dt 



plotted in Figures 4, 8 and 12,. The isollnes of 

~ fd^ j'^(yp-v) are drawn in Figures 5? 9 and 13. The field 
g 

of effective divergence is plotted in Figures 6, 10 and 1H-, 

and isobaric divergence is plotted in Figures 7? 11 and 1?. 

The labeled values of isolines of fd ^ D and 
2 \ .3 g Dt 

7^-yj are times lO”'^ and the labeled values of the 



g 
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fields of divergence are times 10 . In plotting the values 

of Isobaric divergence in Figure 7 the values in four adja- 
cent squares were averaged and then plotted in the common 
point of the four squares.. This was done in an effort to 
smooth the divergence field but was not done in Figures 11 
and 15 because the averaging could mask the effects of short 
waves. 



5.- The Results of the Computations. 

The individual terms on the right hand side of Equation 

-3 

7 have order of magnitude 10 feet per second. A typi- 
cal value of the change of (z^-z) over a twelve hour period 
is -100 feet. When divided by 12 hours or 4.32 x lO^seconds, 
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iz-z) has a value of about t 2.3 x K^T^feet per second. 

“3 

The extreme computed value of (^-z) was 11.6 x 10 feet 
per second. Measuring ^Vp(r-z)\ over a distance of 5 of lati- 
tude each side of the station) gives values of the order 

-5 

of 5.0 X 10 . With a wind of 120 feet per second and an 

* 

angle ^ whose cosine is 0.3 j a typical value of term (b) 

is 1.8 X lo”^ feet per second. The maximum magnitude of 

this term was 20.5 x 10 feet per second. The quantity 

ranges from 11,4 x 10”^ at 30*^ and 60° N to 
S 6 

13.18 X 10” at 45° where it is maximum. When multiplied 

I 

by a typical value of 100 feet per second for v, one obtains 
a value for term (c) of 1.1 x 10“^ to 1.3 x 10”^ feet per 
second. The maximum value of term (c) was I.-87 x 10"^ feet 
per second. Although term (c) is of the same order of magni- 
tude as terms (a) and (b) , it is somewhat smaller in general 
and terms (a) and Cb) give the largest contribution to the 
sum of the three terms of the left side of Equation (7) . 

The value of f^d^ is constant for a given latitude 
circle and ranges from 1,37 x 10 feet at 30° N to 3^66 x 10 
feet at 55° N, while (z-z) ranges from -300 feet to 300 feet. 
Term (d) was positive in all the computations except one,, in 
which it was a small negative value. Of course all compu- 
tations were carried out in the broad scale flow between 
trough and ridge. Term (d) could very well be a negative 
value in a closed high pressure system viiere values of 
(z-z) might be as great as -400 feet. 

Values of isobaric divergence were in general smaller 

8 



than those of effective divergence. However they agree 
fairly well qualitatively.. The magnitude of the effective 
divergence agrees very well with that computed by Landers 
The average absolute value of the effective divergence on 
the three charts was 3.3^) 3.70 and 1.22 all times 10**^ per 
second. A plot of divergence values versus height was made 
by Landers and shows an increase of 0.6 x 10”^ per 

second in the mean absolute value of horizontal divergence 

from 20,000 to 25,000 feet. Extrapolation of his curve to 

-5 

30,000 feet would give a value of 2.7 x 10 per second, and 
the mean absolute value of effective divergence of all the 
computation made in this study is 2.78 x 10”^ per second. 

The values of isobaric divergence have a maximtim magnitude 

-5 

of about 3.5 X 10 per second while the mean absolute 

value is about 1.0 x 10 per second which agrees better 

with the values of divergence as computed by Arnason Ql^ . 

In the following section a statistical approach is made in 

comparing the two fields of divergence and the fields of 

fd£ ^ and -f^ 
g Dt g 

6. Statistical Studies. 




Figures 12 and 13 seem to give the best agreement. To test 
this agreement a correlation coefficient of the two quantities 
was found using values at 35 points in each figure. This 
correlation coefficient r is given by the foibwing equation 

9 



from Hoel C*+J* 



(9) r - coy (xv) 

2y 

Cov(xy) is the covariance of x and y, the stochastic variables, 
and Sjj and Sy are the respective standard deviations, Cov(xy) 
and the standard deviations were computed from the following 



equations 


given in [V3* 


(10) 




(11) 


si 


(12) 


s| - ] 



The mean value of fd^ D was -2,06 x 10^ feet per 

g Dt _3 

second and the mean value of - fd^ ( t7p-y / was 0,1 x 10 

g ^ ^ -3 

feet per second. The standard deviations were 6,38 x 10 
feet per second and the 3*52 x 10"^ feet per second respect- 
ively, The correlation coefficient was found to be 0,31* 
With 90 per cent confidence limits -the critical value of 
r is 0,29 so the correlation coefficient of 0,31 is signif- 
icant at this level of confidence. However, the value of 
r is not large enough to deduce more than that the two 
quantities are related. Some of the reasons for this lack 
of agreement will be discussed later. 

The same procedure as described in the previous para- 
graphs was used to compute a correlation coefficient between 
Figures 10 and 11 which are patterns of effective divergence 
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and isobaric divergence. It was believed that perhaps better 

agreement could be achieved between these two patterns since 

in the process of dividing D !fa by -f some of the errors in 

Dt 

using the geostrophic approximation might be minimized. How- 
ever, the correlation coefficient between Figures 10 and 11 
was only O.IM+ which is significant at 68 per cent confidence 
limits with a sample size of 52 * 

The approximation of -A-. (5-z) by a 12-hour average 

^ t 

would tend to cause a computed effective divergence center 
to be east of its instantaneous position if the future 12- 
hour (z-z) map is used and west of its instantaneous posi- 
tion if the previous 12-hour ("^-z) map is used. The plotted 
values of effective divergence and isobaric divergence were 
examined in an attempt to find some consistent displacement 
which would give a higher correlation coefficient. It was 
believed that the correlation coefficient could be improved 
by shifting the divergence center over Indiana in Figure 10 
into coincidence \-/ith the divergence center over Michigan 
in Figure 11 while the convergence centers over Iowa remained 
in coincidence. Actually, however, there was no improvement 
in the correlation coefficient. A value of 0,07 was com- 
puted for r with a sample size of 31 which is not significant 
at any reasonable level of confidence. 

Sutcliffe has stated that the vortex- tube term in * 
the vorticity equation tends to be offset by Y) 

so that a better approximation to the vorticity equation 
would be : 
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A correlation coefficient \ias computed using tliis 
approximation and the computations from the l^OOZ, 16 
November chart. Using Equations (9)j (10), (11) and (12) 
the correlation coefficient was computed to be -0.204 which 
is significant at 75 per cent confidence limits. The fact 
that the correlation coefficient is negative seems to indi- 
cate that ignoring fC?r y)^ust overcompensate for the vortex- 
tube term. 

Using the equation for isobaric divergence in natural 
coordinates: 



ci4) V, • y 



V Aii. 4- ^ ^ 



it is seen that there are two components of the total di- 
vergence, a longitudinal component along a streamline and 
a lateral component normal to the streamline. It was 
also of interest to determine v/hich of these components 
gives the greatest contribution to the total isobaric di- 
vergence. The rate of change of wind velocity along a 
streamline is most easily measured and was found at all the 
stations at which vorticity computations were made. This 
was done by measuring the change in wind speed along a 
streamline for a distance 5° of latitude across each sta- 
tion. If these values are compared with the effecti-'^e 
divergence, the difference will be the effective lateral 

component. In most cases the lateral and longitudinal 

12 



terms are of opposite sign. In 103 out of 13^3 computa- 
tions (approximately 75/0 the lateral term was the domin- 
ant terra. Therefore, most of the time the sign of the 
effecti’^e divergence is determined by the effective 
lateral component. 

It would be helpful if some relationship could be 
derived between effective divergence and the longitudinal 
component since the longitudinal component can be easily 
measured. The total divergence could then be determined 
simply by measuring ^ . Using Equations (9)> (10), (11) 
and (12) to determine a correlation coefficient between 
and the effective divergence gave a coefficient of 
0.087 vrfiich is not significant at any reasonable level 
of confidence. Therefore, very little about the magnitude 
of the effective divergence can be deduced from the magni- 
tude of the longitudinal divergence. 

7. Sources of Error. 

The discrepancy which exists between computed values of 
2 2 

f^ Djgj. and -fd_f(a.\/) and between effective and iso- 
g Dt g ^ 

baric divergence is due partially to the fact that the vortex- 
tube term and the vertical advection of vorticity have been 
ignored. However, the discrepancy is greater than the values 
to be expected from the terms that have been ignored. Due to 
the lack of actual wind reports at 300 mbs, it vms necessary 
to interpolate between the actual wind observations. Error 
could have been introduced here» Also a 10° error in the 
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observed wind direction and a 5-knot error in ;d.nd speed 
may cause a maxiraura error of 3.0 x 10“^ per second in the 
computed isobaric divergence v/hile the same error in vrf.nd 
direction and speed cause only an error of 0.3^ x 10”^ 

in the effective divergence. However, there is another 
approximation in the use of the vorticity equation. This 
was the determination of ■— Cz-z) as a 12-hour average 
using two consecutive Cz-z) -charts. The magnitude of the 
errors due to this approximation is unkno'tm, 

8, Synoptic Aspects. 

It was difficult to determine the track of the diver- 
gence centers as many of them disappeared or were out of the 
area of computations in the 12-hour intervals^ For example, 
in the 12-hour interval from Figure 6 to Figure 10, three 
divergence centers have appeared where there was previously 
only one. 

The effective divergence patterns give good qualitative 
agreement with the synoptic situation.. Fleagle [2^ has- found 
the 500-mb level to be approximately a level of non-divergence 
and the 300-mb surface is near the level of maximum divergence, 
Landers [5j also found the 300-mb level to be a level of maxi- 
mum divergence and the 500-nb level to be a level of least 
divergence. Arnason's results [^l] did not show such a distri- 
bution of divergence, but it is believed that the isobaric 
divergence as^ computed at 300 mbs may give a good indication 
of rises and falls of the 500-mb surface. In Figure 6 the 
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areas of effective divergence between trough and ridge are 
relatively large and clear cut. With this divergence 
pattern aloft we would expect the trough to continue to 
move eastward, wiiich actually occurs. Figure 7 does not 
gi-^^e this same agreement as there is a convergence center 
centered over station 53^ which is behind the ridge. This 
would indicate retrogression of the ridge which does not 
occur. It may also be noted in Figure 6 that as the strong 
winds east of the ridge at 84° W enter the cyclonically 
curved isohyp ses of wealter gradient there is a change from 
convergence to divergence. On the left hand side of Figure 
6 there appears to be a perturbation leaving the long-wave 
trough as indicated by the divergence center over Colorado, 
As yet there is no actual short-wave trough indicated at 
500 mbs. On the next 500-mb chart, 0300Z 16 November, a 
short-wave trough has appeared east of the long-wave trough. 
Figures 10 and 11 indicate closed divergence centers ahead 
of the short-wave trough and closed convergence centers 
behind it, V/ith the short-wave trough in the broad scale 
flow,, the pattern of divergence is much more complicated and 
the gradients are much greater. On the West Coast in Figure 
10 and 11 there is a center of convergence located ahead of 
a short-wave trough, T4/elve hours later the trough is prac- 
tically non-existant. The divergence center located over 
Oregon may be associated \d.th the trough that appears off 

I 

the West Coast on the 1500Z 16 November chart. In Figures 
l4 and 15 the short-wave trough of Figures 10 and 11 lias 



moved into the ridge and weakened. Consequently the divergence 
and convergence centers that were associated with it are 
not .as clear cut. The patterns have become elongated 
along the contours with considerable convergence east of the 
trough. Therefore, the trough may be expected to decrease in 
amplitude and the flow to become more zonal, which is verified 
on the 03 OOZ, 17 November 500-mb chart* 

Due to the fact liihat the surface maps are two and a 
half hours earlier than the upper-air charts, it is diffi- 
cult to compare divergence patterns with the surface maps. 

\ 

However,, in the few cases studied the low-pressure centers 
at the surface move tovrard the centers of di-'^ergence at 
300 mbs* It is also interesting to note that areas of di- 
‘vergence tend to overlay the fronts at the surface,. 

9 , S\jinmary and Conclusions. 

This investigation of the. isobaric divergence patterns 
at 300 mbs has shown that the technique of computing diver- 
gence from the rate of change of absolute vorticity gives • 
very good agreement in magnitude with the computations of 
isobaric divergence made by Landers [^5^ • computations 

made in this investigation indicate that a mean absolute 
value of divergence at 300 mbs is about 2,8 x 10”^ per second* 
This value of effective divergence acting over a column of 
air 12,000 feet thick for a period of 12 hours will cause a 
fall in the 500-mb surface of about 200.^ feet. The values of 
effective divergence plotted in Figures 6, 10 and 1^- give 

good indications of future areas of rise and fall at the 

16 



500-mb surface. 



The computations of Isobarlc divergence using Equation 
(8) gave a smaller mean value, v;hich gives good agreement 
with the magnitudes computed by Arnason {^ij . Except in 
Figure 7> these values of divergence also give good indi- 
cations of future developments on the 500 -mb surface. 

There were rather poor statistical correlation co- 
efficients between computed values -of fd^ D 

2 \ g 

• V and also between effective and isobaric diver- 

g ■ ~ 

gence,. This poor agreement is believed due largely to ^-errors 
in, and insufficient, wind reports at 300 rabs» Also the 
12-hour approximation of Introduces error. 

The lateral component of the effecti-'/e divergence v/as 
the dominant component in 75^ of the computations made.- 
There was found to be no relationship between the effective 
divergence and either the lateral or longitudinal component 
of dl-'^er gence.. 
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